Impurity and spin effects on the magneto-spectroscopy 
of a THz-modulated nanostructure 
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We present a grid-free DFT model appropriate to explore the time evolution of electronic states 
in a semiconductor nanostructure. The model can be used to investigate both the linear and the 
nonlinear response of the system to an external short-time perturbation in the THz regime. We 
use the model to study the effects of impurities on the magneto-spectroscopy of a two-dimensional 
electron gas in a nanostructure excited by an intense THz radiation. We do observe a reduction 
in the binding energy of the impurity with increasing excitation strength, and at a finite magnetic 
field we find a slow onset of collective spin-oscillations that can be made to vanish with a stronger 
excitation. 



I. INTRODUCTION 

The interest in the properties of radiation induced 
excited states in bulk semiconductors^*^ or a two- 
dimensional electron gas (2DEG))^ and semiconductor 
nanostructurcs^ has grown in the last years parallel to 
advances in high precision experimental techniques and 
sample making skills. Increased brilliance of coherent far- 
infrared (THz) synchrotron radiation^ and sensitive pho- 
toconductivity measurements techniquesi*^ will certainly 
intensify research activities on electronic systems driven 
far from equilibrium and time-dependent phenomena in 
nanostructures. 

We model and investigate here the time evolution of 
a short nanoscale quantum wire or an elongated quan- 
tum dot with one embedded impurity after the system 
with several electrons has been irradiated by a short in- 
tense THz pulse. Corresponding atomic systems have 
been investigated)^ but the present model of an artificial 
two-dimensional (2D) molecule in a GaAs structure al- 
lows us to consider the effects of an external magnetic 
field i? < 2 T, that is commonly available in the labs 
and has large effects on the 2DEG in a AlGaAs/GaAs 
heterostructure, but very small effects on usual atomic 
systems. 

Having seen that an intensive THz pulse can cause non- 
linear effects in quantum dots2i and rings^ and possibly 
change the binding energy of impurities in bulk material? 
we study its effects in the present system on the binding 
energy of the impurity in presence of other interacting 
electrons, and ask ourselves if it is possible to excite col- 
lective spin modes in the system in an external magnetic 
field. THz radiation with the equivalent energy of few 
meV is simply on an appropriate scale with the energy 
spectrum of our system determined by its geometry (fi- 
nite size) and the host material of GaAs. 

Modeling time-dependent phenomena in nanosystems 
in a magnetic field driven out of equilibrium does not 
leave many options to a researcher. The Coulomb in- 
teraction between the electrons is usually of paramount 



interest together with the spin structure of the system. 
This can be emphasized by its finite size even though the 
effective g-factor may be small as is the case in GaAs. 
The well established time-dependent Hartree-Fock ap- 
proximation can be adequate in many cases^^ but one 
would hope for more accurate spin structure delivered 
by a time-dependent local spin-density approximation 
(LSDA).^ It has to be kept in mind that both the inclu- 
sion of the magnetic field in a LSDA^" and its extensions 
to dynamic systemsii are nontrivial tasks and are not 
fully completed to date. In our approach we shall use 
the ground state parameterization of the exchange and 
correlation functionals and ignore the small corrections 
due to current terms in the functionals in a magnetic 
fieldiS 

Commonly, both in calculations for excited atoms and 
moleculesSiiiii or nanostructures j^iS a spatial grid has 
been employed. In this work we resort to recently devel- 
oped grid-free implementation of LSDA^^^iiSiiLiS in the 
hope to be able to describe complex dynamical systems 
for a longer period of time. The grid-free methods, as will 
become clear below, result in the transformation of the 
relevant equations to compact matrix equations of not 
too large size. On the other hand the matrix elements 
are usually much more complicated than in grid methods, 
even though most of them are known analytically. This 
compact matrix structure of the problem lends itself very 
well to the matrix and vector handling of modern com- 
puter languages and facilitates parallelization that does 
neither require lare bandwidth nor high communication 
speed between the nodes of a CPU cluster. 



II. MODEL 

A. The ground state 

To calculate the time evolution of the system requires 
the knowledge of the initial state, its ground state, at 
a particular time t = tg, which within the LSDA is de- 
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scribed by the Kohn-Sham equations 

H\a) - (i/o + H, + V^ + i/i„t) \a) = £„|a), (1) 

where is the Hamihonian for the 2DEG in the per- 
pendicular magnetic field B — Bz, that is confined by 
the parabolic potential Vconf('') = m*ujQr^ /2. A more 
general shape of the system is attained through the ad- 
ditional confinement potential breaking the rotational 
symmetry of Hq and modifying the radial confinement 
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^ ap cos {pip) + Vo exp {-jr'^), (2) 



where Pmax is an even integer to ensure confinement. The 
eigenstates \a) of the Fock-Darwinii^ Hamiltonian Hq are 
used as a mathematical basis in which the eigenstates \a) 
of H are expanded, as a result Eq. can be transformed 
into a matrix eigenvalue equation. A suitable truncation 
of the basis {|a)}, that can be varied to improve the 
accuracy, then leads to a numerically tractable eigenvalue 
problem. 

In this LSDA mean-field type of an approximation 
the effective interaction Hmt is composed of the direct 
Coulomb repulsion potential, or the Hartree term 



^H(r) 



r — r' 



7dr\ 



(3) 



and the potential describing the effects of the exchange 
and correlation 
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(nea;c[f^T,"i,^])|n„=n„(r), (4) 



where na-{r) is one of the two spin (a =|, I) components 
of the total electron density n = n-\ + ni. The Zeeman 
energy of the electrons is included in with the appro- 
priate effective g-factor g* and the Bohr magneton, fiB, 
= ±(l/2)g>BB. 
In a large 2DEG the natural length and energy scales 
are the magnetic length / = y^hc/{eB) and the cyclotron 
energy hujc — heB/{m*c), and instead of the electron 
density it is convenient to define the local filling factor 
i/(r) = 2TTpn{r). In the confined system described by Hq 
the appropriate length scale is a = l^Joj^/Tl, an "effec- 
tive magnetic length" , where Vl — u)"^ -\- Auiq is the new 
characteristic frequency of the system. We shall thus use 
a modified filling factor appropriate for a finite system 



^{v) = 27ra^n(r), 



(5) 



and the polarization C(r) = ["•t('") ~ n^w 
variables to express the exchange and correlation poten- 
tials in. This will allow us to have a unified descrip- 
tion of the system in a finite or vanishing magnetic field. 

In these new variables the exchange and correlation 
potentials are^ 



d d 
VxcA = -^{vexc) + (1 - O-^^xc 
ov oQ 

d d 
VxcA = -^iye^c) - (1 + O-^^xc- 



(6) 



In the magnetic field the exchange and correlation energy 
density is interpolated between the infinite and the no- 
field limits as^"'^ 



^xci^^ C) = + el^{V, 0(1 - (7) 

with f{D) = {■iv/2) + lv'^ and e;^(i>) = -0.782\/^eV(fcO- 
The low-field energy functional is interpolated between 
the spin-polarized and the unpolarizcd limitsS^ 



4(^,0 = exc(i>,0) + r(C) [e.c(i>,l) 



with r(c) 



(i + C)^/^ + (i-C ) 

23/2 _ 2 



=(^,0)] 

3/2 _ 9 



(8) 



The exchange and correlation contributions are then sep- 
arated e,j;c{i',0 — ^x{i^iC) + £c(^', C); aiid the exchange 
contribution for the unpolarizcd system is expressed as 
ea;(i>, 0) = — [4/(37r)]\/^e2/(fcZ), while for the polarized 

one it is e^iv,!) = -[4/(37r)]\/2i>eV(fcO- We use the 
parameterization of Ceperley and Tanatar for the corre- 
lation contribution^-^ 



1 + aix 



+ aix + 022^ + a^x-^ 



Ryd-- 
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where x = y/r^ — {2/i'y^^{a/a*gy/'^, and a*g is the effec- 
tive Bohr radius. The optimized values of the correlation 
parameters have been found by Ceperley and Tanatar 
using a Monte Carlo model of a 2DEG,^- 

So stated, the V^co-'s are complicated functionals of 
the the spin densities, or D and C, and thus a com- 
mon method is to solve the resulting Kohn-Sham equa- 
tions on a spatial grid. We instead follow the alter- 
native example of a number of researchers to use the 
mathematical basis {|a)} to cast the LSDA functionals 
into matrix expressions in order to implement a grid- 
free LSDAiSiiSiiLia Briefly, the method relies on the fol- 
lowing steps: In the end we need the matrix elements 
of the exchange and correlation potentials, {a\Vxc,iT\P), 
so we start by constructing (aji'l/S). The polarization 
C has to be treated as a combination of two functions, 
f — i'l — S'l and g — i/^^ . Due to the completeness of 
the basis we then have (a|C|/3) — X]7('^l/l7)(7lffl/^)' ^'^'^ 
to evaluate the latter matrix elements we need to find 
a unitary matrix U that diagonalizes the matrix i/, i.e. 
V — \S diag(Ai, • • • , A„)U"'", where Ai are the eigenvalues 
of V. All functions of the matrix v are now calculated 
according to 



/p] = Udiag(/(Al),•••,/(A„))U^ 



(10) 



Corresponding methods are used for more complex func- 
tionals of V and C^. 

The matrix elements of the electron density n or the 
effective filling factor v can be written in terms of a gen- 
eralized one-electron overlap integral and the matrix el- 
ements of the density operator p 



^Pqv J dr <j)lir)<j)*p{r)<Pg{r)<j)p{r). 



(11) 
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In the ground state the density matrix is constructed 
from the expansion coefficients of electronic states in 
terms of the basis states, the occupation of the states 
expressed by the Fermi distribution, and the chemical 
potential /i ensuring the conservation of the number of 
electrons, pqp = '^-ffi^-y — t^)CypC-yq- We use here the 
density matrix pqp in anticipation of the situation when 
the system will be subjected to a strong external time- 
dependent perturbation and quantities like the single- 
electron energy spectrum, and the occupation described 
by the Fermi distribution will have no meaning. The an- 
alytic expression for the generalized one-electron overlap 
integral can be found in Appendix. The matrix elements 
of the direct Coulomb interaction and the confinement 
potential have been published elsewhere i^iiS^ 

B. Time-dependent excitation 

At i = <o the Hamiltonian of the system acquires time 
dependence, H{t) — H + W{t), caused by an external 
perturbation of a finite duration 

W{t) = Vtr^^''^ cos {Np(t)) exp {- sr^ ^Tt) 

sin {ui it) sin (ujt) 9 {tt ~ u! it). (12) 

Since the perturbation is of arbitrary strength we have 
to assume the task of solving the equation of motion for 
the density operator 

iMtp{t)^[H + Wit),pit)]. (13) 

The structure of this equation is inconvenient for numer- 
ical evaluation so we resort instead to the time-evolution 
operator T, defined by p{t) = T{t)poT^{t), which has 
simpler equation of motion 

iht{t) = H{t)T{t) (14) 
-inf+{t) = T+{t)H(t). 

In each time step the density operator p changes and thus 
also the DFT Hamiltonian which is a functional of p. 
Similar nonlinear problem is encountered when the time- 
dependent Kohn-Sham equations are solved directly on 
a spatial gridi^ We apply a similar procedure here dis- 
cretizing time and using the Crank-Nicholson algorithm 
for the time integration resulting in the coupled equations 

p{t + Af) = T{At)p{t)T+{At) (15) 
{l + ^i/[p;. + A.]}T(At).{l-f 

with the initial condition, r(0) = 1, already used on the 
right hand side of (|15|l . In each time-step the equations 
have to be iterated until self-consistency is attained. This 
is performed in the truncated Fock-Darwin basis {\a)} for 
the corresponding matrix version of the equations (|15|l . 

The accuracy and the fine-tuning of the model has been 
tested for the ground state with comparison to the ex- 
act results calculated for a circular parabolically confined 



quantum dot)2& and in the time-evolution by convincing 
us that the model can reproduce the Kohn results for 
the absorption of a Far-infrared radiationiSLSS In addi- 
tion, the normalization condition for the density matrix 
is tested regularly during the integration of the time- 
dependent problem. 

C. Specification of model parameters 

In order to model a short GaAs quantum wire with 
or without an embedded impurity and several electrons 
we select the parameters, hojo = 3.37 meV, and assume 
g* = 0.44, and m* — 0.067m. The impurity potential 
is selected to be Gaussian Q with Vb = — 10 meV, and 
a^7 = 1.0. The Gaussian shape is used here in order 
to limit the size of the basis during the CPU-intensive 
time integration of the problem. The wire shape of the 
system is attained by setting a2 = —0.7 and = for all 
j ^ 2. The charge density for 10 electrons in the system is 
shown in Fig. ^ with and without an embedded impurity. 
At _B = T the length of the system is approximately 
200nm, and the difference between the charge densities 
with and without an impurity in the system turns up in 
Fig.ffla as a fairly well localized impurity state. This can 




FIG. 1: Color online. The electron density of the interacting 
system without an impurity (a), the difference in the electron 
density of system with and without an impurity (b), and the 
density for a system with an impurity (c). i? = T, Vo = — 10 
meV, T = 1 K, iV = 10, and a « 13 nm. 

be confirmed by the Kohn-Sham energy spectra displayed 
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in Fig. 121 where the main change due to the addition of 
the impurity potential into the system is the lowering 
of the energy of the highest almost degenerate doubly 
occupied state by approximately 1 meV. The calculations 




FIG. 2: Color online. The Kohn-Sham single-electron energy 
spectra for the interacting electron system with (□) and with- 
out (*) an impurity for B = T (a), and B = 1 T (b). The 
integer a is a quantum number labeling the states. The lower 
horizontal line (yellow) represents the chemical potential ^ for 
system with an impurity, and the higher one (blue) indicates 
^ for a system without an impurity. T = 1 K and A*' = 10. 

reported here are for i? = T, and i? = 1 T. In order 
to avoid numerical problems with the determination of 
the eigenvectors for exactly degenerate eigenvalues in the 
diagonalization of the Hamiltonian we actually use B = 
10"^ T instead of B = T. 

For the radiation pulse ifT^ of duration just over 3 
ps^ we select F = 2 ps~^. The envelope frequency cor- 
responds to fiuji = 0.658 meV, and the base frequency 
is equivalent to fiLO — 2.63 meV. Meanwhile, we select a 
|iVp| = 1 and s = in ((T^ to represent a dipole radia- 
tion. But, due to the finite cut-off in the basis \a) the 
matrix elements of the excitation can never be purely of 
the dipole type. We choose the orientation of the initial 
radiation pulse (|12|l to be either parallel, ||, to the long 
axis of the system or perpendicular to it, _L. The finite 
asymmetric cut-off of the basis can introduce a tiny mix- 
ing of the perpendicular and the parallel excitation at 
vanishing magnetic field. In fact, the basis cut-off is per- 
formed by specifying the highest radial quantum number 
riniax used in the calculation, by doing so an asymmetry 
is introduced between the number of negative and pos- 



itive angular quantum numbers M . This asymmetry is 
made unimportant with respect to the energy spectrum 
and the wave functions of the system by selecting a high 
enough nmax, but the representation of a dipole poten- 
tial is not quite so clean. A simple alternative method is 
not easy since the effective asymmetry does depend on 
the magnetic field. 

We integrate the equations for the time-evolution oper- 
ator H15|) over a time interval of 50 ps using an increment 
At = 0.0025 ps and evaluate the expectation values for 
the coordinates x and y for each spin direction, (x)]-, etc. 
We can thus identify oscillations perpendicular or parallel 
to the long axis of the system, talking about perpendic- 
ular or parallel detection. Resonances in the collective 
oscillations of the electron spin densities are sought by 
Fourier transformation of the time series in the interval 
5 — 50 ps, when the external excitation has been turned 
off and switch-on effects have vanished and the system 
has reached a "steady state" . 



III. RESULTS 

A. The B = T case 

In Fig. 131 the Fourier power spectra of the time series 
for the coordinates x and y are used to represent the ex- 
citation spectra of the system without an impurity. The 
II polarized excitation is seen to introduce a tiny _L os- 
cillation (1/10 the amplitude) as has been explained ear- 
lier, but at the same frequency. This simple excitation 
spectrum does not change much with increasing excita- 
tion strength Vt- A second harmonic shows up in the 
parallel oscillation, but more interestingly the height of 
the main peak does not grow linearly with the excitation 
strength when Vta > 1 meV. The excitation then be- 
comes nonlinear and as V^a = 5 meV we see that it shifts 
slightly to lower energy. Before commenting on the red 
shift of the peak with increasing Vt we show what hap- 
pens in the system with an embedded impurity in Fig. 
0] The main peak in the longitudinal oscillation (Fig. 
EJj) is slightly blue shifted compared to the system with 
no impurity. The reason is readily explained by a glance 
at the Kohn-Sham energy spectra in Fig. 12 the pres- 
ence of the impurity potential enlarges the energy gap 
around the chemical potential /i. Instead of a spatially 
fairly extended state just below /i there is now a localized 
state at lowered energy. Again, with increasing excita- 
tion strength Vt we see a red shift of the resonance in 
the parallel oscillations. At strong excitation we can not 
refer back to the Kohn-Sham spectra of the ground state, 
but we can expect that the energy pumped into the sys- 
tem makes the states just below the gap more extended, 
thus reducing the effective energy gap. The perpendic- 
ular oscillation in Fig. is more complex than before. 
Again, there is a resonance at the same energy as in the 
parallel mode that is slightly red shifted with increased 
excitation, but there now appear resonances at higher 
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FIG. 3; Color online. The Fourier power spectrum of the 
time series for the expectation values of the coordinates per- 
pendicular (_L) (a), and parallel (|j) (b) to the long axis of the 
system. The excitation is parallel to the same axis. Division 
factors for the stronger excitations are indicated in the fig- 
ure. No impurity is in the system. B = T, T = 1 K, and 
N = 10. 



FIG. 4: Color online. The Fourier power spectrum of the 
time series for the expectation values of the coordinates per- 
pendicular (_L) (a), and parallel (||) (b) to the long axis of the 
system. The excitation is parallel to the same axis. Division 
factors for the stronger excitations are indicated in the figure. 
Impurity is present in the system. B = T, T — 1 K, and 
iV = 10. 



energies indicating the presence of more localized modes 
in the system, an extra excitation mode located at the 
impurity site. This is again supported by a look at the 
Kohn-Sham energy spectra in Fig.[21i, at the changes just 
above the gap caused by the impurity. 

The case of the initial excitation pulse aligned perpen- 
dicular with the long axis is review in Fig.[5l We see that 
again in the _L detection we have the peak structure en- 
countered in Fig.^ as could be expected, but the sHght 
presence of a || excitation in the initial pulse is not of 
the dipole kind so mostly higher order modes are excited 
in the || direction. Exactly this, is further displayed in 
Fig. IH^, where, in addition, the relative strength of the 
response is also visible. In contrast, the strong confine- 
ment in the perpendicular direction of the system or the 
wire ensures the _L response to be not strongly dependent 
on the main excitation polarization, see Fig. Eb. 

To clearly demonstrate the collective oscillations in the 
system we show in Fig.Elthe induced electron density for 
the system with and without an embedded impurity for 



the two polarizations of the excitation pulse. The two 
figures Fig. [7^ and b for the _L excitation show that the 
ends of the system oscillate in the _L direction while the 
oscillation around the impurity is at least at this mo- 
ment parallel to the long axis. In case of no impurity 
in the system the oscillation of the center part is more 
perpendicular. When the excitation is applied || to the 
system (see Figs. [7|; and d) the main oscillation mode 
is also along the system of the fundamental dipole type 
as can be seen by comparing the induced density to the 
total density in Fig. ^ In this case the oscillation seems 
to be more extended to the whole system when there is 
an impurity present in it, but when there is no impurity 
present the center of the system displays a tiny bit of a _L 
mode superimposed on its motion. In fact, the induced 
density can show higher order oscillation modes that we 
can not discern with our time series of the expectations 
values for the the coordinates x and y, which are most 
sensitive to dipole oscillations occuring in the system^*S 
The red shift of the collective resonance for the system 
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FIG. 5: Color online. The Fourier power spectrum of the 
time series for the expectation values of the coordinates per- 
pendicular (_L) (a), and parallel (|[) (b) to the long axis of the 
system. The excitation is perpendicular to the same axis. 
Division factors for the stronger excitations are indicated in 
the figure. Impurity is present in the system. B — T, 
T = 1 K, and iV = 10. 



FIG. 6: Color online. Comparison of the Fourier power spec- 
trum of the time series for the expectation values of the coor- 
dinates parallel (||) (a), and perpendicular {!.) (b) to the long 
axis of the system in the case of an excitation perpendicu- 
lar or parallel to the same axis. A multiplication factor is 
indicated in the figure. Impurity is present in the system. 
Vta = 1.0 meV, B = T, T = 1 K, and iV = 10. 



with an embedded impurity is further explored in Fig. |S| 
in comparison with the response of the system without 
an impurity. Clearly, the resonance is red shifted and 
broadened with increasing excitation strength Vt pointing 
towards a reduction in the binding energy of the impurity. 

Comparison between the sub figures of Fig. |S1 reveals 
that with increasing excitation strength more energy 
flows into a simple dipole oscillation of the system when 
an impurity is present in it. This is confirmed by the 
induced electron density in this case showing large oscil- 
lations of the electron density around the impurity site 
parallel to the system, and mostly in harmony with the 
oscillations in the rest of the system, reflected in a single 
resonance. The broadening of the resonance is indicative 
of a finer structure that we are unable to resolve in our 
calculations. 

The case is different with the perpendicular oscillation 
shown in Fig.|51 where in the end a strong excitation will 
lead to many modes becoming active. So, even though 



the parallel collective oscillation remains simple the per- 
pendicular one is generally more complex. 



B. The B = 1 T case 

The collective oscillations in the finite magnetic field 
_B = 1 T are for short times after the initial excitation 
pulse very similar to the oscillations encountered in the 
system in no magnetic field as could be guessed by com- 
paring the Kohn-Sham energy spectra for the two cases, 
see Fig. |21 It has though to be kept in mind that the 
Lorentz force now ensures a quick communication be- 
tween the parallel and the perpendicular modes, they 
are not independent. But more happens: The finite size 
of the system together with the small Zeeman energy, 
and the finite temperature T = 1 K establish a tiny spin 
polarization in the ground state, i.e. the two spin den- 
sities do not have exactly the same spatial dependence. 




FIG. 7: Color online. The induced density at time t — 50 
ps for ± excitation for a system without (a), and with an 
impurity (b), and for || excitation without (c), and with an 
impurity (d). The excitation strength Vta — 1.0 meV (a), 5.0 
meV (b) 10-2 meV (c), 1.0 meV (d). B = T, T = 1 K, and 
N = 10. 



The spatially varying exchange force in the LSDA gives 
a different restoring force to the spin densities. All these 
factors lead to a slow flow of energy into the individual 
oscillations of the spin densities as can be seen in Fig. ^| 
This flow is so slow that after 50 ps the system has still 
not entered a "steady state" as in the case of i? = T. 
We are still observing switch-on effects for the collective 
oscillations of the spin densities after 50 ps. In the spin 
oscillation we see present both high and low frequency 
components but the lack of a steady state prevents us 
from making a Fourier analysis. Just to avoid any mis- 
conception we should mention that the time increment 
of At = 0.0025 ps gives us 20,000 points in the range of 
50 ps and at a finer level of resolution the spin oscillation 
looks very smooth with no apparent noise like might be 
hastily inferred by looking at Fig. [TUl 

Even more amazing is the fact that to begin with this 
spin oscillation increases with the strength of the exci- 
tation Vt but, as is seen in Fig. ^| for a still stronger 
excitation it again decreases. Stronger excitation of the 
system prevents the energy in flowing into the collective 
spin oscillations, but directs it to the center-of-mass mo- 
tion and the other collective charge oscillations. This 
behavior is definitely a nonlinear phenomena that can 
not be observed in a model of linear response. 

The presence of the impurity also quite effectively di- 
rects more energy into the spin oscillation as is seen in 
Fig. 1101 This can be explained by a stronger spatial vari- 
ation of the exchange force around the impurity. 
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FIG. 8; Color online. The Fourier power spectrum of the time 
series for the expectation values of the coordinate parallel (||) 
to the long axis of the system with or without an impurity in 
the case of an excitation along the same axis. B = T, T = 1 
K, and N = 10. 



The regularity of the collective spin oscillations can be 
seen in Fig. II II where simply both spin densities are seen 
almost overlapping at i = 50 ps. The peak structure 
of the spin densities closely follows the structure of the 
ground state charge density displayed in Fig. ^ 



8 



1 1 Excitation 


No imp. 

Imp. 


J_ Detection 




Vta=10"^meV 




(a) 





3 4 

E (meV) 



1 



3 4 5 

E (meV) 





No imp. 




Imp. 


1 1 Excitation 


J_ Detection 




Vta=1 .0 meV 




(b) ^ 











No imp. 

Imp. 


1 1 Excitation 




± Detection 




Vta=5.0 meV 






(c) ^ Jv/v^JWf 


,1 





012345678 

E (meV) 



FIG. 9; Color online. The Fourier power spectrum of the 
time series for the expectation values of the coordinate per- 
pendicular {!.) to the long axis of the system with or without 
an impurity in the case of an excitation along the same axis. 
B = T, r = 1 K, and TV = 10. 
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IV. DISCUSSION AND SUMMARY 

Without the external magnetic field present the collec- 
tive charge oscillations excited in the system by the short 
THz pulse do reach a "steady state" condition within 1 
or 2 picoseconds after the excitation is switched off al- 
lowing a Fourier analyzes to discern the modes present. 



FIG. 11: Color online. The two components of the spin den- 
sity at the time point t = 50 ps in a system with an impurity. 
The spin densities are very similar in shape and are superim- 
posed here on one another with different shades. Vta — 10~^ 
meV, B = 1 T, r = 1 K, and = 10. 
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At low excitation level we observe a linear regime, but 
at higher excitation nonlinear effects show up. The ex- 
citation gap in the system is reduced indicating a less 
localized single-electron states and in particular the ef- 
fective binding energy of the impurity is reduced. 

At finite magnetic field (i? = 1 T) a slowly increasing 
collective oscillation of the spin densities has not reached 
a "steady state" condition after 50 ps. As the spin os- 
cillation is directly coupled with the collective charge os- 
cillations a Fourier analyses is not possible in this case 
for either component. The excitation of the collective 
spin oscillations is much more effective when an impu- 
rity is present in the system, since then the initial spin 
densities a.t t — have a larger spatial variation. The 
more amazing fact is that with increased intensity of the 
initial excitation the onset of the spin oscillation can be 
effectively shut down. We contribute this to the flow of 
the energy nonadiabatially pumped into the system be- 
ing more easily directed to the oscillations of the charge 
density in this case. 

For the results shown in this paper we have always used 
the same type of excitation pulse. What would happen if 
we use a longer pulse with several oscillations or of a dif- 
ferent form. We can not claim that we have thoroughly 
investigated this question, but our limited experiment- 
ing with this seems to indicate that the pulse form is 
indeed important for the heights of the resonance peaks, 
but does not affect their location strongly. We are pump- 
ing energy into the system by the excitation in varying 
amounts depending on the peak shape and duration and 
it can be expected that this also affects the initial distri- 
bution of the energy between the modes. By removing 
the system further from equilibrium than we do we are 
sure that stronger nonlinear effects can occur with con- 
siderable change of peak locations also. We already see 
the precursors of this in the redshift of the peaks con- 
nected to the impurity embedded in the system. 

We are aware that our choice of an impurity potential 
represented by a Gaussian well does not correspond to 
unscreened Coulomb impurities that might be situated 
in the plane of the 2DEG in a nanosystems. This form 
has only been selected in order to keep the tnmcation 
of the Fock-Darwin basis within reasonable bounds in 
order to enable the "long-time" integration of the time 
evolution operator to extend to 50 ps with 20,000 steps, 
within each one the self-consistency was attained with 
further 4 loops to properly update the density operator. 
With these requirements in mind it should be reiterated 
that the present model was not designed to describe any 
known experimental device, but rather it was constructed 
to give qualitative insight into time-dependent phenom- 
ena that can take place in a nanosystem that is irradi- 
ated by short and intensive THz pulse, in response to the 



preparations underway in several experimental groups to 
study these phenomena. We only study the behavior of 
the system for a short time after the application of the ex- 
citation and neglect all possible damping and dissipation 
effects that would be caused by the coupling of the sys- 
tem to its immediate environment. These are certainly 
very important effects that await further investigation. 
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THE GENERALIZED ONE-ELECTRON 
OVERLAP INTEGRAL 

In order to reduce the memory size needed to store the 

generalized overlap integral it is convenient to split it into 
two parts by applying the completeness of the basis once 

■' 5 

(16) 

where an analytic expression for the wave functions corre- 
sponding to the Fock-Darwin basis states \a) = \Ma, ria) 
gives 

^ (^ly.+v+i. /n^ + \M^\\ (np + \M, 

la—Olp—Ols—O ^ ^ 

n5+\Ms\\ / |M^| + |M^| + |M^| 
ns-is J \ 2 

2\ [i\Mc\+\Mp\+\Ms\)/2]+ia,+ip+is+l 



with 



+ ia+ip + 'i'5 ! 



rial 



2(|M<,|+i)/2'W 7r(|M„|+na)' 



(17) 
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